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ABSTRACT 

Spitzer Space Telescope infrared measurements are presented for 24 members of the TW Hydrae association (TWA) . 
High signal-to-noise 24 /xm photometry is presented for all of these stars, including 20 stars that were not detected by 
1 IRAS. Among these 20 stars, only a single object, TWA 7, shows excess emission at 24 /*m and at the level of only 
40% above the star's photosphere. TWA 7 also exhibits a strong 70 /im excess that is a factor of 40 brighter than the 
■ stellar photosphere at this wavelength. At 70 fj,m, an excess of similar magnitude is detected for TWA 13, though no 
' 24 fim excess was detected for this binary. For the 18 stars that failed to show measurable IR excesses, the sensitivity 
| of the current 70 /im observations does not rule out substantial cool excesses at levels 10-40 x above their stellar 
— . ■ continua. Measurements of two T Tauri stars, TW Hya and Hen 6-300, confirm that their spectacular IR spectral 
' energy distributions (SEDs) do not turn over even by 160 /im, consistent with the expectation for their active accretion 
disks. In contrast, the Spitzer data for the luminous planetary debris systems in the TWA, HD 98800B and HR 4796A, 
q ■ are consistent with single-temperature blackbody SEDs and agree with previous IR, sub-millimeter, and millimeter 
5h ' measurements. The major new result of this study is the dramatic bimodal distribution found for the association in 
c/3 , the form of excess emission at a wavelength of 24 /im, indicating negligible amounts of warm (> 100 K) dust and 
C$ • debris around 20 of 24 stars in this group of very young stars. This bimodal distribution is especially striking given 
^ ' that the four stars in the association with strong IR excesses are > 100 x brighter at 24 /xm than their photospheres. 
Clearly, two terrestrial planetary systems, HD 98800B and HR 4796 A, exist in some form. In addition, there are at 
least two active accreting objects, TW Hya and Hen 6-300, that may still be forming planetesimals. The remaining 
stars may possess significant amounts of cold dust, as in TWA 7 and TWA 13, that have yet to be found. 
Subject headings: infrared: stars — circumstellar matter — planetary systems: formation 



1. INTRODUCTION 

While calibrating the Infrar ed Astronomical Satellite 
(IRAS) lAumann et all JlQSl and ICilletti (fl98^ dis- 
covered large far-infrared excesses around four nearby 
stars: Vega, Fomalhaut, (3 Pictoris, and e Eridani. Re- 
lease of the IRAS Point Source Catalog led investiga- 
tors to many other young, nearby dwarf stars display- 
ing excess emission at m i d- an d far-IR wavelength s (e.g 



Walker fc WolstencroftJ 119881 Ide la Reza et all 
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1992; Mannings & Barlow 



1989; 



1998; 



Zuckerman fc Becklinl 11993). As a result, numerous 



investigations of these objects from the ground have 
added critic al insights concerning cir cumstellar ma- 
terial (e.g.. Uavawardhana et all Il998t iKc 
1998; IKenvon fc Bromlevi 120021 EOT 



voerner 



The Hubble 



et al.l 
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Space Telesco pe (HST) and the Infrared Space Obser- 
vatory (ISO) fflabing et alJl200Tt TSpangler et aHl20uTl 
IMever fc Beckwithl I2000T and references therein) have 
also given us key insights into specific objects, including 
images taken in the near-IR showing scattered light from 
thin disks of dust that are shown to abs orb starlight and 
re-em it at much longer wavelengths (e.g- JSchneider et alJ 
1999). Thus, in the last 20 years a body of observa- 
tions, measurements, analysis and theory have accumu- 
lated that leave no doubt about the reality of terrestrial 
planetary material around young stars. Now, investiga- 
tions of planetary debris systems (PDSs) are underway 
with t he new tools provid ed by the Spitzer Space Tele- 
scope ijWerner et al.ll200^) . 

We have initiated a program using Spitzer to explore 
a unique and varied sample of late-type stars that are 
nearby, ~40-100 pc, and show one or more indications 
of young age. Our total sample of over 120 stars has been 
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prioritized into three groups: (1) the first 19 stellar sys- 
tems identified by common space m otion as members of 
the TW Hydrae a ssociation (TWA; iKastner et aT]ll997t 
iWebb et "all fl99 9"l . (2) ~60 stars that show significant 
lithium absorption indicating young ages, noting that all 
of the TWA systems easily meet this requirement and 
therefore overlap this sample, (3) a group of ~ 60 stars 
with less certain ages based on either detection of X-rays 
or Ha in emission. 

A combination of special circumstances dictate 
that we start our exploration with the TW Hy- 
drae association. First, the estimated age of the 
association ranges from 8—10 million years (Myr) 
and the ages of these stars are essentia lly equa l 
ijStauffer. Hartmann. fc Barrado v Navascuesl [l995). 
Second, the TWA members are among the closest young 
stars available. Third, most association members lie in 
excellent, low-background regions of the sky for far-IR 
observations. IRAS has already shown that four of the 
brightest IR excesses known are found in the TWA. 
Additional important factors are that five of the TWA 
stars have Hipparcos distances averaging around 50 pc, 
their spectral types range from AOV to M5V with most 
objects being of M or late K spectral type, and almost 
all members show signs of T Tauri activity. TW Hydrae 
itself has been classified as a cl assic T Tauri s ystem 
despite its apparent iso lation ijHerbig. G. H.I Il978t 
iRucinski fc Krautterlll9 83). In contrast, and within the 
same association, the brightest and most dramatic PDS, 
HD 98800, shows no evidence for current accretion of 
primordial material onto the star and has served as the 
most direct proof of ongoing terrestr ial planet formation 
Ce.g.. lLow. Hines. fc SchneiderllT999j) . 

2. OBSERVATIONS 

Observations of 24 members of the TW Hydrae asso- 
ciation were made using Spitzer. These include 14 single 
stars and binaries that are unresolved at 24 /im and five 
binary systems at least partially resolved by Spitzer at 
24 /im. The Multi-band I maging Photometer for Spitzer 
(MIPS; Rieke et al.H2004|) provided measurements or up- 
per limits of the flux densities (F v ) at 24 /im and 70 /im 
for the entire sample. MIPS 160 /im photometry was 
also obtained for the four IR-excess systems discovered 
by IRAS (i.e., TW Hya, Hen 3-600, HD 98800B, and 
HR 4796 A). The observations were made between 2004 
January 30 and June 21, with each target observed in all 
of the intended wavelength bands during a single visit 
by Spitzer. Table 1 lists the stars in the association ob- 
served with MIPS along with basic data from previous 
measurements. Six other system s have been ident i fied as 
memb ers of the TWA (see e.g., ISong et al.ll2003t IWebbl 
l2000a|) since the final planning of the Spitzer /MIPS ob- 
servations, but were not included in the target list be- 
cause of time limitations. 

The MIPS data are reduced using the data analysis 
tool ( DAT) of the MIPS instrument team ijGordon et all 
120051) . The images are corrected for distortion, regis- 
tered, and mosaicked to compensate for the movements 
of the spacecraft and cryogenic scan mirror mechanism 
(CSMM) within MIPS that are employed to produce a 
series of spatially dithered images on the three detector 
arrays. Aperture photometry is then performed on the 
calibrated, mosaicked images. Table 2 summarizes the 



MIPS bandpasses, aperture photometry parameters, and 
the calibration factors adopted for the data reduction. 
The data have been color corrected using the relevant 
tables in the Spitzer Observer's Manual. 

The 24 /im data were collected using the standard 
MIPS point source observing template that obtains 14 
images of the target at various standard locations on the 
128 x 128-pixel Si: As detector array during each cycle of 
the template. Either 3 s or 10 s integration times and 
1-4 complete template cycles were used to ensure high 
signal-to-noise ratio (S/N) measurements of the targets. 
Total observation times were estimated for each target 
so that a S/N>10 measurement was achieved at the ex- 
pected 24 /xm flux level of the star's photosphere, and all 
stars were detected to at least this S/N. Indeed, the me- 
dian S/N is ^160, though the internal photometric cal- 
ibration is estima ted to be 1-2% in this bandpass (e.g., 
IRieke et al.ll200l . 

For the aperture photometry of the 24 /im data, an 
aperture 15" in radius (r) was typically used along with 
a sky annulus centered on the aperture and having an 
inner radius of ~63" and width of 25". The large sky 
annulus was chosen to avoid as much flux as reasonable 
from the Spitzer plus MIPS point spread function (PSF) 
of the target. A larger aperture (~37"5 radius) was used 
for four binary TWA systems that are at least partially 
resolved at 24 /im (TWA 8, 9, 13, and 15). The rela- 
tive brightnesses of individual components of the bina- 
ries were estimated by comparing the fluxes measured 
in r = 2"5 apertures centered on the components. For 
TWA 9 and 13, a third source of 24 /im flux fell within the 
large photometric apertures. The flux from these sources 
were subtracted and are not included in the values for F v 
at 24 /im. 

Only HD 98800 was bright enough that the core of the 
PSF saturated even in the minimum observation time 
allowed at 24 /im. In this case, the flux in the satu- 
rated core was estimated by measuring the flux in the 
wings of the PSF. This was done for the entire sample 
as well. The exercise on the non-saturated stars showed 
that the 24 /im flux in the saturated portion of the PSF 
for HD 98800 could be recovered to within 5%. After 
correcting for the saturated core, the 24 /im flux density 
was measured for HD 98800 in the usual manner. 

Observations, data reduction, and aperture photome- 
try for the 70 /im data are similar to those at 24 /im. 
However, the pixel scale of the 70 /im wide-field optical 
train (~10" pixel -1 ) dictates that none of the binaries in 
the TWA can be resolved except for TWA 19. Therefore, 
we have used a uniform photometry aperture (3-pixel ra- 
dius) and sky annulus (radius = 4-6 pixels) for the stars 
observed in the 70 /im wide-field mode. The chosen pho- 
tometry aperture is very large compared to the < 1" 
pointing accuracy of the telescope, eliminating the am- 
biguity in both the target position on the array and the 
determination of flux density upper limits for undetected 
sources. 

TW Hya, HD 98800, and HR 4796A were observed in 
the 70 /im narrow-field mode to prevent saturation by 
these bright IR sources as this optical train halves the 
pixel scale. For these observations, an r = 5"67 aperture 
was chosen. Instead of an annulus centered on the target 
to estimate the sky contribution within the photometric 
aperture, rectangular regions on either side of the target 
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were utilized because of the narrow, 16-pixel width of the 
usable portion of the detector array. 

Because of the degraded sensitivity and higher instru- 
mental noise of the 70 /im Ge:Ga detector array relative 
to pre-launch estimates, it was not practical to obtain 
measurements accurate enough to detect objects at the 
expected levels of the TWA photospheres. As a result, 
the 70 /im observations were tailored to detect excesses 
at flux levels >10x the expected 70 //m photosphcric 
flux density with a reasonably high degree of confidence 
(>3-4cr). Unfortunately, most of the TWA observations 
occurred before the bias level for the 70 /im detector 
was reduced. During the period before the bias change, 
the overall noise in the detector was greater and there 
were more pronounced changes seen in background lev- 
els across array modules. Because of these instrumental 
effects, a conservative criterion of 5a above background 
was adopted for a detection at 70 /tm. 

All four targets observed at 160 /zm were detected at a 
S/N ranging from ^6-160 using the standard MIPS pho- 
tometry observing template. In all cases, an aperture of 
r = 64" was used for the photometry. The sky contri- 
bution within the photometric aperture was estimated 
using rectangular regions on either side of the target. As 
with the 24 /im and 70 /im data, the 160 /mi photome- 
try is corrected for the loss in flux caused by the use of 
apertures of finite size not encompassing the total flux 
from the PSF. The multiplicative aperture corrections 
to a hypothetical aperture of infinite extent are listed in 
Table 2. 

The 160 /im data were also corrected for the blue fil- 
ter leak discovered in this bandpass after the launch of 
Spitzer. This leak allows light at ~2 /im in wavelength 
to reach the detector which is sensitive to near-IR ra- 
diation. Tests made during flight have shown that the 
contaminating flux is about 15 x the flux received from a 
star's photosphere at 160 /im assuming a Rayleigh- Jeans 
law between the near and far-IR. In all four cases in 
the TWA where 160 /im measurements were made, the 
expected photospheric 160 /im flux density is <Cl% of 
the measured signal. Given the bright far-IR excesses of 
these stars, the filter leak corrections are small. 

In addition to broad-band photometry, observations 
of TW Hya and HD 98800 were made using MIPS 
in its SED mode. These observations were made on 
JD = 2453395, roughly a year after the MIPS photom- 
etry of TW Hya was obtained and about seven months 
after the photometry of HD 98800. In this observing 
mode, the CSMM brings light from the telescope to a 
reflective diffraction grating that also performs the func- 
tions of a slit and collimator. The collimated, dispersed 
light is then focused onto the 70 /tm detector array. The 
result is a very low-resolution (R — A/AA ~ 15-25) first- 
order spectrum covering a wavelength range of 52-97 /im. 
The slit/grating has a width of two pixels (~ 20") and 
the dispersion is 1.7 /im pixel -1 . An unusable portion of 
the 70 /im array restricts the slit length to ~ 2.'7. Also, 
a dead 4 x 8-pixel readout at one end of the slit further 
reduces the effective length of the slit where the full spec- 
tral range can be sampled. 

Observations in SED mode involve chopping the 
CSMM between the target and a region 1-3' away to 
sample the background, and using small spacecraft moves 
to dither the target between two positions along the slit. 



Data reduction is again handled by the DAT in essen- 
tially the same manner as for 70 /im broad-band imag- 
ing. The result is co-added "on-" and "off-target" im- 
ages. After the background image has been subtracted 
from the target, the spectrum is extracted using an ex- 
traction aperture of set width (8 pixels). Flux calibration 
of the TW Hya and HD 98800 spectra was accomplished 
using a MIPS SED-mode observation of a Boo made on 
JD = 2453196. The 50-100 /an continuum of a Boo is 
assumed to be cx A -2 and F v — 14.7 Jy at the effective 
wavelength of the 70 /im filter bandpass. The observa- 
tions and spectral extractions of a Boo and the two TWA 
systems were identical and the same DAT calibration files 
were used for every target. 

3. RESULTS 

Table 3 summarizes the MIPS photometric results for 
the TWA. Several stars in the association are known 
optical variables and we give the epoch of the start 
of the observations in the second column of the ta- 
ble. The 24 /im, 70 /im, and 160 /im flux densities are 
given in units of milli-Janskys (mjy; where 1 mjy = 
10~ 26 erg cm~ 2 s _1 Hz -1 ). Uncertainties in the flux 
densities are the square root of the quadratic sum of the 
measurement and calibration uncertainties, noting that 
most relative uncertainties are much lower. Also listed 
in Table 3 is the S/N of the detections in each MIPS 
band. Although some of the S/N values are extremely 
high, the uncertainties quoted for the photometry in all 
three bandpasses are dominated by uncertainties in the 
absolute flux calibration of the instrument. These are 
estimated to be ~10% for the 24 /im band and ~20% 
for the 70 /im and 160 /im bandpasses. Upper limits as- 
signed in Table 3 to non-detections at 70 /tm are at the 
level of 3a over the noise of the background. 

Our primary goal for this observational program 
is to detect and measure IR excesses in the TWA. 
This requires that good estimates of F v be made for 
the TWA stellar photospheres at the effective wave- 
lengths of the MIPS bandpasses. For all of stars 
in the sample, stellar photospheres were derived us- 
ing a grid of Kurucz l | K urucz| | 1979D and "Next Gen" 
l)Hauschildt. Allard. fc Baronll99 9'l photospheric models 
and fitting available optical photometry combined with 
near-IR flux densities compiled from the Two-Micron All 
Sky Survey (2MASS) Point Source Catalog. Tempera- 
ture is a free parameter in the fitting, but in all cases solar 
metallicity and logg = 4.5 are assumed. The NextGen 
models are known to yield better fits to late-type main 
sequence stars than Kurucz models, and were used for 
all TWA members except for HR 4796A and TWA 19A, 
which have much higher effective temperatures than the 
rest of the sample. The best-fitting model temperatures 
(T») of the TWA are listed in Table 1. 

3.1. 24 /iw Excesses 

Comparison of the 24 /im data for the TWA with the 
2MASS -ftT s -band flux densities of these stars shows that 
the association can essentially be divided into two pop- 
ulations. Figure 1 plots the 24-to-2.17 /tm flux density 
ratio against the flux density in the -fCj-band. In all cases, 
the 2MASS photometry is consistent with the stellar pho- 
tospheres emitting the near-IR light. The four stars with 
IR excesses discovered by IRAS are about a factor of 100 



4 



Low et al. 



brighter at 24 /im relative to K s than the other stars in 
the sample. In fact, for the stars with flux ratios <0.02, 
the data are consistent with the 24 //m flux being solely 
from the photosphere. As a guide, a dashed line in Fig- 
ure 1 marks the limit of F(24 /j,id)/Fk s if both band- 
passes fall within the Rayleigh- Jeans regime of the stel- 
lar continuum. Since the TWA members are typically 
spectral type K and M, the K s band falls close to the 
emission peak of the photosphere. As a result, the stars 
typically lie above the Rayleigh- Jeans limit in Figure 1. 

Figure 2 further explores the empirical relation be- 
tween the 24 /im flux density with the near-IR. The stars 
with previously known IR excesses or T» > 4500 K are 
omitted from Figure 2, and the flux ratio shown in Fig- 
ure 1 is now plotted against the temperature of the model 
stellar photosphere. The majority of the 19 stars plot- 
ted form a locus that is offset from the F(24 ^,vo)/Fk b 
ratios expected from blackbodies. Spectra derived from 
the NextGen models generally give a good approximation 
to the observed flux ratios, although the TWA members 
have not reached the main sequence. The unweighted av- 
erage flux ratio (F(24/xm)/ 'F Ks ) = 0.015 with an RMS of 
only 0.004. The concentration of the points close to the 
expected values from the NextGen models in Figure 2 
suggests that we are seeing no IR excess at 24 /im for 
most of these late-type, young stars. The MIPS photom- 
etry is consistent with the 24 /im flux emanating from the 
photosphere for most of the members of the TWA. One 
example is TWA 6, where t he MIPS 24 /im mea surement 
supports the contention bv lUchida et al.l l|2004|) that the 
5-20 /im spectrum of this st ar taken with Spitze r and the 
Infrared Spectrograph (JRS: lHouck et afll2004|) is indeed 
the spectrum of the photosphere. 

Although the NextGen models yield a much better 
brightness estimate of the photosphere at 24 /im than 
blackbodies in this temperature range, Figure 2 shows 
that the flux ratio is possibly more dependent on the 
model temperature than is actually observed. Note, 
TWA 12, 13B, and 17 fall well below the NextGen curve. 
These deviations and the intrinsic scatter of the data 
are intriguing, but the limited number of objects in the 
sample complicate a detailed investigation of their pos- 
sible causes, especially considering that only tempera- 
ture was varied in constructing the models and that the 
TWA stars are younger than the nominal zero-age main 
sequence stars described by the NextGen models. In the 
case of TWA 17, the 24 /im flux density is measured to 
be only 1.5 ± 0.2 mJy and the source is detected at a 
S/N — 13 with MIPS. Therefore, we do not confirm the 
possible excess at 12 and 18 /im reported from ground- 
based observations by IWeinberger et afl (|2004fl . 

The measurement of F(24 /im)/fjf s for TWA 7 sug- 
gests that it possesses a 24 /im excess that contributes 
~40% to the total brightness of this system at this wave- 
length. Except for the four IRAS detections, TWA 7 
has a higher flux ratio in Figures 1 and 2 than all of the 
other TWA members. The next highest ratio belongs to 
TWA 9B, but its deviation from the expected flux of the 
photosphere is less than half of that shown by TWA 7, 
and is no larger than the flux deficits at 24 /im mea- 
sured for TWA 12, 13B, and 17. Since the 2MASS and 
MIPS measurements are not simultaneous, an alternate 
explanation for the apparent IR excess in TWA 7 could 
be the variability of the star which is not uncommon for 



T Tauri stars, although there are no published data sug- 
gesting that this star is variable. Also, the clustering of 
most measurements around (F(24 /im)/Ff S - s ) is a clear 
demonstration that most stars in the TWA do not, in 
fact, exhibit significant variability. 

3.2. 70 fim Excesses 

Regardless of the strength of a possible 24 /im excess 
in the observed continuum of TWA 7, this object is de- 
tected at 70 /im at a flux level ~40x brighter than the 
photospheric emission estimated for the star at this wave- 
length. Figure 3 presents the ratio between the observed 
and predicted 70 /im flux densities. Only six of the 20 
systems (TWA 19 is resolvable at 70 /im) are detected 
with MIPS at 70 /im. TWA 13 joins TWA 7 and the 
four IRAS sources in having a detected 70 /im excess. 
TWA 13 cannot be resolved with Spitzer at 70 /im, so it 
is unknown which of the two components produces the 
excess, or if both components of the binary contribute to 
the observed flux. 

Figure 3 shows the same basic division in the TWA 
sample as is found at 24 /im. The 70 /im excesses for 
TWA 7 and 13 are >10x fainter than the other detec- 
tions. One-sigma upper limits are also shown in Figure 3 
for the remainder of the sample. As discussed in §2, the 
sensitivity limits of the 70 /im array render detection of 
the photospheres of the TWA sample impractical. Ex- 
cesses similar to those of TWA 7 and 13 cannot be ruled 
out for many members of the association, and it is clear 
from the data that it is not necessary for an object to 
show an excess at 24 /im to find one at 70 /im. For in- 
stance, there is no hint from the 24 /im measurement 
of TWA 13 of an excess at longer wavelengths. There 
are three stars where a 70 /im excess at >10x the flux 
density of the photosphere can be ruled out at the 3cr 
confidence level: TWA 2, 5A, and 8A. Therefore, the 
minimum range in strength for 70 /im excesses in the 
TWA is greater than a factor of 300. 

3.3. Spectral Energy Distributions 

The SEDs of the six systems in the TWA with con- 
firmed IR excesses are shown in Figure 4. It is readily 
apparent that this small sample of young stars exhibits a 
rich variety of IR properties. The T Tauri stars TW Hya 
and Hen 3-600 show no hint of their SEDs turning over 
even out to 160 /im. In fact, the MIPS 160 /im mea- 
surement for TW Hya reveals that the spectrum of the 
presumed optically thick accretion disk continues to rise 
into the infrared. Similar to TW Hya, the MIPS data 
for Hen 3-600 are consistent with the earlier IRAS re- 
sults and show that the 70 and 160 /im flux densities are 
roughly equivalent. 

In contrast, HD 98800 A and B do not show the ex- 
treme activity associated with T Tauri stars. How- 
ever, HD 98800B generates a huge IR excess that 
is wel l fit by a single blackbody with T = 160 - 
170 K i|Low. Hines. fc SchneiderlT99a limes et all2004|) . 
whereas component A lacks an easily measured excess 
even though the total luminosities of the two compan- 
ions are virtually identical. The MIPS photometry is 
consistent with earlier results. Likewise, the IR excess 
of HR 4796A produces an IR SED that can be fit by a 
single blackbody, but it is si gnificantly cooler than found 
for HD 98800B (T = 110 K; lJuraet al~I IT993). With the 
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addition of the Spitzer observations, we fit this excess 
with a T = 108 K blackbody. Given that HR 4796A 
is an AOV star, it is not surprising that it has an IR 
excess because of the finding that younger A stars have 
a higher probability of possessing a 24 pm excess and 
havin g stronger therma l excess emission than older A 
stars ijRieke et all 120051. It is of note that HR 4796 A 
has by far (by ~10x) the largest ratio of 24 /im excess 
emission-to-photo spheric emission of the 265 A and B 
stars examined by iRieke et alJ l|2005|1 . 

Figure 5 presents the MIPS SED-mode observations 
of TW Hya and HD 98800B. These low-resolution spec- 
tra do not show any spectral features of high equivalent 
width in cither object. Arcturus was used to calibrate 
the spectra and they are consistent with the 70 /im flux 
densities derived from the photometry. The IR spectral 
slopes of TW Hya and HD 98800B are quite different, 
with the continuum of TW Hya being much redder than 
that of HD 98800B. A power-law having a spectral index, 
a, where F v oc u a , is adequate to describe these data. 
For TW Hya, a = -0.33 ± 0.04, and is consistent with 
the finding that the 160 /jm flux density is higher than 
the 70 /im measurement. The slope of a single power- 
law fit to the continuum of HD 98800B between 52 and 
97 /im yields a = 1.12 ± 0.04. The blackbody fit with 
T BB = 160 K to the PDS of HD 98800B is also shown 
in Figure 5, and its slope within this spectral region is in 
good agreement with the SED-mode data. 

3.4. PDS Luminosity and Mass 

Our primary reason for calculating the stellar and IR 
excess luminosity is to measure the fractional luminosi- 
ties of those stars that show direct evidence of PDS for- 
mation and to provide well defined upper limits for those 
stars that fail to show an excess. Table 4 lists the stellar 
luminosity, L*, derived by numerical integration, and the 
fractional luminosity, Lm/ X*, where Lm is the luminos- 
ity of the infrared excess that is the result of numerically 
integrating the flux density in excess of the predicted 
level of the photosphere. Direct measurements of the 
distance to TWA members are limited to only five sys- 
tems (Table 1). We have assigned a distance of 55 pc, the 
distance corresponding to the unweighted average of the 
parallaxes measured by Hipparcos excluding TWA 19, to 
the other stars in the sample. In general, L* is higher 
than expected for the TWA given their spectral types 
because these young, rapidly evolving stars have higher 
luminosities, lower temperatures, and larger radii than 
when they reach the main sequence. 

Conservative upper limits for Lm/ X* have been as- 
signed to objects listed in Table 4 that show no evidence 
for a 24 /im excess by assuming that the la upper lim- 
its estimated for the 70 pm flux densities represent the 
brightness of a possible IR excess plus photosphere at 
this wavelength. It is also assumed that any possible IR 
excess is well approximated by a blackbody with tem- 
perature, Tpjjs- The upper limit on the temperature 
of a possible cold PDS is given in Table 4 and is set so 
that Tpds an d the assumed value for F(70 /im) raise 
F(24 /im) by no more than 10%; the adopted uncer- 
tainty in the absolute flux calibration of the MIPS mea- 
surements at 24 /im. These upper limits on temperature 
do not preclude small amounts of warmer dust within 
systems that are still consistent with the 24 pm mea- 



surements, but they do constrain the properties of cir- 
cumstellar material able to produce substantial excesses 
at 70 /xm and be undetected by Spitzer. 

The mass of circumstellar dust responsible for the lu- 
minosity of IR excesses observed in the TWA, or that 
could produce the upper limits found for Tpp>s and 
Ltr/L*, c an be deduced followin g the inferences of 
Uura et alJ 1(1991 and lChen fc Jural p00l|l . That is, the 
minimum mass in dust can be estimated by 

M d > (16/3)Tr(L m /L*)pR PDS (a) 

l(Jura et alJ 119951) . where (a) is the average radius of a 
dust grain, Rpds is the minimum distance from the star 
where the grains are in radiative equilibrium, and p is the 
mass density of the grains. The equation above assumes 
that the dust is radiating from a thin shell at a distance 
Rpds from the illuminating star and that the debris disk 
is optically thin. Also, it is assumed that the dust grains 
are spherical and have cross sections equal to their geo- 
metric cross sections. In Table 4, we tabulate Md using 
the values for averag e grain size and density adopted by 
IChen fc Jural po¥ ((a) = 2.8 pm; p = 2.5 g cm" 3 ). 
The minimum dust mass responsible for detected IR ex- 
cesses range from ~ 10 25 g for HR 4796A to ~ 10 23 g 
for TWA 13. For HD 98800B, the assumption that the 
debris disk is optically thin is not valid, accounting for 
the much smaller mass es t imate than that determined by 
lLow. Hines. fc Schneiderl l(1999(K 

The mass estimates in Table 4 for the TWA are conser- 
vative lower bounds for the total circumstellar dust mass 
in these systems since the MIPS observations are not sen- 
sitive to colder material located further from a star and 
emission from larger particles. For example, there is evi- 
dence for material around TWA 7 colder than the ~80 K 
dust detected by Spitzer. Observations at 850 pm by 
iWebbl l)2000aH detected TWA 7 with F v ~ 15 mJy This 
flux density is about a factor of seven higher than the 
extrapolation to 850 pm of the ~80 K blackbody d eter- 
mined by th e 24 and 70 pm excesses. IWebbl l(2000a(l and 
I Webb et alJ (|2000bfl interpret the sub-millimeter emis- 
sion as coming from 20 K dust with grain sizes on the 
order of a few hundred microns. As a consequence, the 
mass estimate for the dust disk around TWA 7 derived 
by the sub-millimeter measurement is ~ 10 3 x greater 
than that estimated from the IR excess. 

4. DISCUSSION 
4.1. Lnfrared Properties of the TWA 

The main result of the Spitzer-MIPS observations of 
the TWA is to show that most stars in the association 
exhibit no evidence for circumstellar dust and that the 
24 pm data place severe limits on the presence of warm 
(T > 100 K) dust within these systems. Although the 
TWA includes four remarkable IR-excess systems, there 
are no other objects that come close to their 24 pm out- 
put relative to the observed near-IR photospheric emis- 
sion. The bimodal distribution of warm d ust in the 
TWA, identified bv iWeinberger et all 1(20041) at shorter 
wavelengths, is confirmed at 24 pm for a larger sample 
of objects (see Figure 1). 

The bimodal nature of 70 pm excesses is not as pro- 
nounced as at shorter wavelengths, but there is still a gap 
of about a factor of 10 between the strengths of the 70 pm 
excess of Hen 3-600 and the Spitzer detections of TWA 7 
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and 13. There is also a range of at least ~ 300 between 
the brightest and faintest of the 70 /im excesses in the 
TWA relative to their photospheres at this wavelength. 
The luminosity of the observed IR excesses range from 
~ 25% (TW Hya and HD 98800B) to - 0.1% (TWA 13) 
of L*. 

It is difficult to fully characterize the IR excesses of 
TWA 7 and 13 given that they only begin to appear 
at ~24 /im. In the case of TWA 13, there is also the 
complicating factor that the system is a binary and the 
distribution of the material responsible for the excess is 
unknown. However, calculating a color temperature for 
the two objects yields similar results: T c ~ 60-80 K. 
Assuming that the excesses are well fit by blackbodies, 
both systems have a fractional IR luminosity of ~1% of 
that measured for HD 98800B; L IR /L* = 2.0 x 10~ 3 and 
1.7 x 10~ 3 , for TWA 7 and 13, respectively. Although 
the lack of observed excess 24 /im emission constrains 
the temperatures of possible PDSs for most of the 70 /im 
non-detections to be < 100 K, the MIPS observational 
upper limits cannot rule out IR excesses with fractional 
luminosity in the range of 10 -3 -10 -4 . Even at the es- 
timated young age of the association it may already be 
more likely to detect cooler material, further away from 
the stars, than warmer material that results in an excess 
at 24 /um or shorter wavelengths. Clearly, at 70 /im wave- 
length, large amounts of additional Spitzer time will be 
needed to detect and measure the cool dust and debris 
around these stars. 

The variety of the IR properties in the TWA is remark- 
able given the assumption that its members have similar 
ages. At 8-10 Myr after the formation of the TWA, there 
are still two systems (TW Hya and Hen 3-600) accret- 
ing presumably primordial dust and they have enough 
circumstellar material to reprocess a large fraction (10— 
30%) of the stellar luminosity into the infrared. Also in 
the association are two multiple-star systems, HD 98800 
and HR 4796, that have a prominent PDS around only 
one of the stellar components and the debris is relatively 
warm (T PDS = 100-170 K). TWA 7 shows evidence for 
a weak 24 /im excess with the detected material being 
cooler (Tpds ~ 80 K) and located in orbits > 7 AU 
from the star. The remaining 19 stars, including TWA 7 
and 13, have apparently been cleansed to a high degree 
of dust with T > 100 K. This variety of properties sug- 
gests that dust is swept from regions within several AU 
(see Table 4) of the star very early in the system's evolu- 
tion. The fact that stars in the association show no sig- 
nature of warm dust while other members exhibit strong 
T Tauri activity implies that the transition period from 
the T Tauri stage to the more quiescent state observed 
for most of the TWA K- and M-type stars, is very short. 

An efficient mechanism for the removal of dust is 
thought to be the formation of planets that sweep up 
dust and larger debris in the regions that they orbit. 
iWeinberger et all l)2004|) suggest that planetary forma- 
tion was rapidly completed in the TWA at least in the 
region where terrestrial planets would be expected to 
form. This led to the rapid disappearance of dust in the 
inner several AU of these systems, leaving only the stellar 
photosphere to be detected at wavelengths of 24 /im and 
shorter. Strong stellar winds during the T Tauri stage 
and the Poynting-Robertson drag are also mechanisms 
that can decrease the IR emission from a system by rid- 



ding the environment of small dust grains. Of course, in 
most TWA systems, and within 10 Myr of their forma- 
tion, dust destruction mechanisms must more than com- 
pensate for the dust production mechanism of collisions 
between larger bodies to clear out the terrestrial planet 
region. HD 98800B and HR 4796A are likely cases where 
their debris systems indicate that collisions between bod- 
ies are still important and produce dust close enough to 
the stars to result in large excesses at 24 /zm. 

Nearly all members of the observed sample are 
X-ray emitters, supporting the assertion that stel- 
lar winds in these young, low-mass stars may be 
an important dust destruction mec hanism (see e.g., 
iHollenbach. Yorke. k. Johnstone! |2000|) . We list the X- 
ray-to-stellar luminosity ratios, Lx/L*, for TWA stars 
in the last column of Table 4. The X-ray data are drawn 
from the ROSAT All Sky Survey C atalog and Lx is cal - 
culated in the manner described bv lSterzik etaD |l999). 
iChen et alJ l)2005fl find an apparent correlation between 
IR and X-ray luminosity for their sample of F- and G- 
type stars with ages spanning 5-20 Myr in the Scorpius- 
Centaurus OB Association. This possible correlation is 
in the sense of brighter 24 /im excesses tendi ng to be as- 
sociat ed with fainter stars in X-ray flux, and lChen et alJ 
( 2005) suggest that the stronger stellar wind implied by 
the chromospheric activity may help explain the general 
lack of young systems with 24 /im excesses. Given that 
only three stars in the TWA sample have confirmed PDSs 
that produce measurable excesses at 24 /im, it is not pos- 
sible to test the validity of the X-ray-IR flux correlation. 
However, HD 98800B, HR 4796A, and TWA 7, are all un- 
derluminous in X-rays compared to the average observed 
X-ray emission for the TWA. 

4.2. Comparison of the TWA with the Young Solar 
System 

We can compare the Spitzer observations of young 
stars with models fo r the early So lar System. F ollow- 
ing the discussions in lGaidoa l)1999(l and lJural l)2004|) . we 
assume that the rate of dust production directly scales as 
the rate of lunar cratering during the Late Heavy Bom- 
bardment. The rate of dust production, A'Idust, is given 
by: 

Mdust = Mo (l + 0e*) , 

where tg is the look-back time in 10 9 yr (Gyr), Mq is 
the current rate of dust production and (3 and r are 
fitting constants su ch that (3 = 1.6xl0~ 10 and r = 
0.144 ljChvballT99lj) . In this model, the dust produc- 
tion rate has been approximately constant during the 
past 3.3 Gyr, but was as high as ~ 10 4 the current value 
when t 9 = 4.6. We set M = 3xl0 6 g s" 1 based on 
the z odiacal light in the Solar System ijFixsen & Dweld 
2002), implying that the early Solar System may have 
had a dust production rate of 3xl0 10 gs" 1 . In a model 
where this quantity of dust is produced far from the star 
and then loses angular momentum under the action of 
the Poynti ng-Robertson effect, it is straightforward to 
show fe.g.. lGaidoslll999i) that the luminosity of the dust, 
L u , is given by the expression: 

T MdustC 2 

L v — . 

V 
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For the TWA, young M-type stars would have lumi- 
nosities at 24 /im of 2xl0 18 erg s _1 Hz -1 and, for dis- 
tances of ~50 pc, we expect F l/ (24 fj,m) k 700 mJy, 
vastly greater than what is observed. This result suggests 
the absence of terrestrial planet-forming environments 
around these stars. However, in M-type stars where the 
luminosity is relatively low and the stellar wind rate is 
relatively large, dust grains mainly lose angular momen- 
tum by stellar wind drag rather than through Poy nting- 
Robertson drag (see lJural2004UPlavchan et alJl2005|) . In 
this case, and if M W i n d denotes the stellar wind loss rate, 
then: 

j Mdust 

M wind v 

Although we do not know M W i n d for young M-type stars, 
we can extrapolate from their X-ray emission since the 
winds are likely to be driven by the same hot corona 
which produces the X-rays. It is plausible that young 
M-type stars have mass loss rates 10 2 greater than the 
current Solar wind loss rate of 2xl0 12 g s _1 , and we es- 
timate that M W ind — 2xl0 14 g s _1 . For main sequence 
early M-type stars, we may adopt L* = 0.05L Q . There- 
fore, we expect that at 24 /im, L v = 2xl0 15 erg s _1 Hz -1 
and a predicted flux at 50 pc of only ~0.7 mJy. A 24 /im 
excess this small is beyond the sensitivity of our current 
measurements, and it is possible that planet forming ac- 
tivity may still be occurring around the young M-type 
stars in the TWA even though we do not detect an in- 
frared excess. 

The TWA gives us an important example of the rapid 
evolution of the circumstellar material around generally 
low-mass stars with ages of ~ 10 Myr that is in line 
with time scales of terrestrial planet formation inferred 
for our solar sys t em from rad i ochem istry of meteorites 
ijYin et all 12003 iKleine et all 12002]) . The diverse IR 
properties of the association also suggest that the evo- 
lutionary processes relevant to planetary formation and 
accretion disk dispersal may occur at different rates even 
for stars with similar spectral types. For most of the sys- 
tems observed, the lack of emission from dust within a 
few AU of the stars leads to two possible conclusions. Ei- 
ther, (1) the early conditions around pre-main sequence 
M and late K stars prevent the formation of terrestrial 
planets in most cases, or that (2) terrestrial planet build- 
ing has already progressed to the point where the mass 
is in the form of planetesimals. 
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TABLE 1 
The TWA Sample 



Star a 


Other Name 


Spectral Type 


Ref. b 


distance 

(PC) 


V 


Ref. d 




m f 

2V 
(K) 


TWA 1 


TW Hva 


K8e 


\ 


56 i 7 


11.07 


1 9 


7.30 


3700 


TWA 2AB 


CD-29°8887 


M2e+M2 


2 




11.42 


1-4,6,8,9 


6.71 


3500 


TWA 3AB 


Hen 3-600 


M3e+M3.5 


2 




12.04 


1,2,6,7 


6.77 


3300 


TWA iAR 


STLLJ yoouu 


K4+K5 


□ 
o 


11 ZL U 


o.oy 


1 ( , _LU li) 






TWA 5 A 

_L VVii Jn 




M3e 






11.72 


1 C Q 1 f, 

J. o ,o, y , .LU 


6.75 


3400 


TWA 6 




K7 


4 




12.00 


1 _4 fi io if; 17 


8.04 


4000 


TWA 7 




M 1 


4 




r l.UU 


1 —A Q 1 fi 

i — ii ±,Ujy , J-U 


on 


oouu 


TWA SA 

J- VVii OA 




M2e 


4 




13.30 


1—3 6 12 


7.43 


3300 


TWA 8B 




M5 


1 








9.01 


3100 


TWA 9A 


CD-36°7429A 


K5 


4 


50 ± 6 g 


11.32 


1-6,9,12 


7.85 


4200 


TWA QR 


pn r K(\°7A c >QVK 
Kj VJ — OU I 4:ZyD 


M 1 


4 


Kfl -L- fig 


14 10 




q it; 


'^nn 

OJUU 


TWA in 

-L V V -TV -L \J 




M2 5 


4 




12.70 


1 2 4 6 12 


8.19 


3400 


TWA 1 1 A 


J.LLV, '-t ( 


AO 


5 


67 i 2 


5.78 


9— ZL fi_Q 1 9 


5.77 


9250 


TWA 12 




M2 


6 




13.60 


2,4,12,16 


8.06 


3000 


TWA 13A 


CD-34°7390A 


Mle 


1.6 




12.10 


1,4,12,16 


7.50 


3300 


TWA 13B 


CD-34°7390B 


Mle 


1.6 




12.40 




7.46 


3200 


TWA 14 




MO 


7 




13.80 


2-4,12,16,18 


8.50 


3300 


TWA 15A 




Ml. 5 


7 




14.10 


2,4,16,18 


9.68 


3500 


TWA 15B 




M2 


7 




14.00 




9.56 


3500 


TWA 16 




M1.5 


7 




12.30 


2,4,12,16,18 


8.10 


3700 


TWA 17 




K5 


7 




12.70 


2,4,5,12,16,18 


8.98 


4000 


TWA 18 




M0.5 


7 




12.90 


2,4,5,12,16,18 


8.84 


3800 


TWA 19A 


HD 102458A 


G5 


7 


104 ± 15 


9.10 


2-5,7,12,16 


7.51 


6000 


TWA 19B 


HD 102458B 


K7 


7 


104 ± 15 


11.90 


2,4,12,16,18 


8.28 


4000 



a AB designates binaries that are unresolved by Spitzer at A = 24 fim. 

b References for spectral types a nd i dentifications as me mber s of the TWA — fl)]R i ueim^i &: Krautten 11983|)j (2)[de la Rs za et alj 119891) : 

(3) IGregorio-Hetem etaTl 119921) : C41 IWebb et~al1 119991): (5) IHouU Jl983) : (6) ISterzik et all 1199lf 777) IZuckerman et alJ 120011) . 
c Data are from the Hipparcos Catalog IPerrvman et al.H1997l) . 

d Sources a nd references for pr evious optical -IR photometry — (T l FTbrres et all 120001): ( 2) iReidl 120031) : (3) Hi pparcos-Tvcho IESAI1997I) : 

(4) 2MA SS ICutri et aljl200a) : (5) DENI S lEpchtein et al1[T99a) : (6)ljav awardhan a et all 11 999); (7) IRAS [Bcic hman et al JI1988I) ~f8') 
UCA C1 IZacharias|l200q) : (9) Yale-SPM fPlatais et alJH998B: (10) | Low. Hines, fc Schneiderl 119991): (11) UCA C2 <ZachariasH20dlh (12 ) 
GSC ILasker et alj|19901); (13) ISoderblom et alj 119981); (14) IPrato et alj 120011) ; f 15) IGehrz et aTT 119991) : (16) IWeinberger et all fco04l) ; 
(17) USNO A2.0 IMonettll9*9Sf l: (18) IZuckerman etloT<200ll) . 

c Apparent K a magnitude is from the 2MASS Poi nt Source Catalog, except for HD 98800, where the i^-band magnitudes for both 
components are from lLow. Hines. fc Schneiderl y.999). 

Te mperature of the best-fit Kurucz or Nex t Gen model to the stellar photometry. 
g Scc Mamaick, Lawson. & Fcigclson (2000) for a discussion of the Hipparcos parallax measurement for TWA 9. 
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TABLE 2 
MIPS Bandpasses 



Band 




AA 


scale 




Phot. Apt. 


A b 




W 


(H 


(" pixel- 1 ) 




(radius in ") 




24 fxm 


23.7 


4.7 


2.55 


1.042 x 10~ 3 


15 or 37.5 


1.15 or 1.06 


70 /jm wide field 


71 


19 


9.85 


14.9 


29.6 


1.35 


70 (im narrow field 


71 


19 


5.24 


52.7 


29.7 


1.35 


160 (im 


156 


35 


16 


1.0 


64 


1.33 



a Conversion factor from the calibrated instrumental units of the DAT output to mjy. 

b Factor to correct aperture photometry to that expected from a hypothetical aperture of infinite extent. 



TABLE 3 
MIPS Observations 



Star 


Epoch 


24 (im 


(S/N) 24 M m 


70 jim* 


(S/N) 70 Mm 


160 fim 


(S/N)i 60 Mm 




(JD-2453000.0) 


(mjy) 


(mjy) 


(mjy) 




QC OQO 


9970 -1- 930 


> 1U 




1 I U 


fie: 70 -L 1 Q1 ft 
DO / u nz 10±U 


1 fin 


TWA 9AR 


35.384 


OH 1 -I- 9 n 


200 


< 40 








Hen 3-600 


38.268 


1650 ± 170 


> 10 4 


700 ± 140 


140 


740 ± 190 


6 


HD 98800AB 


177.665 


8500 ± 350 


1000 


6260 ± 1250 


200 


2100 ± 420 


67 


TWA 5A 


38.255 


21.1 ± 2.1 


220 


< 9 








TWA 6 


134.825 


5.7 ± 0.6 


240 


< 8 








TWA 7 


35.398 


30.2 ± 3.0 


260 


85 ± 17 


22 






TWA 8 


35.905 


14.5 ± 1.5 


160 


< 11 








A 




11.1 ± 1.5 












B 




3.4 ± 0.4 












TWA 9 


37.870 


9.3 ± 1.1 


160 


< 34 








A 




6.0 ± 0.6 












B 




3.3 ± 0.3 












TWA 10 


55.726 


5.1 ± 0.5 


190 


< 8 








HR 4796A 


35.928 


3030 ± 303 


> 10 4 


5160 ± 1100 


240 


1800 ± 360 


55 


TWA 12 


177.907 


5.7 ± 0.6 


50 


< 29 








TWA 13 


177.680 


17.6 ± 1.8 


150 


27.6 ± 5.9 


13 






A 




8.8 ± 1.3 












B 




8.8 ± 1.3 












TWA 14 


55.677 


4.1 ± 0.4 


27 


< 25 








TWA 15 


55.701 


2.9 ± 0.3 


72 


< 12 








A 




1.4 ± 0.2 












B 




1.5 ± 0.2 












TWA 16 


55.714 


6.6 ± 0.6 


45 


< 18 








TWA 17 


55.743 


1.5 ± 0.2 


13 


< 13 








TWA 18 


55.762 


2.3 ± 0.3 


20 


< 9 








TWA 19A 


55.689 


10.4 ± 1.1 


100 


< 27 








TWA 19B 


55.689 


4.6 ± 0.5 


44 


< 27 









a Upper limits listed are the detected signal within the photometric aperture plus 3a. If the flux measured within the aperture 
is less than the background level measured for the aperture, the upper limit listed is 3a. 
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TABLE 4 

Properties of Dust Disks in the TWA 



Star 


Tpds (K) 


Rpds* (AU) 


M d (10 23 g) b 




Lir/L, 


Lx/U 


TW Hya 








0.24 


2.7 x 10" 1 


2.7 x 10" 3 


1 WA 2AB 


< 85 


> 7.0 


0.27 


0.39 


< 2.1 X 10 


5.9 x 10 -4 


Hen 3-600 d 








0.35 


1.2 x 10 _1 


6.1 x 10 -4 


HD 98800B d 


160 


2.2 


> 29 c 


0.53 


2.2 x 10 _1 


7.1 x 10~ 4 


TWA 5A 


< 145 


> 2.3 


0.01 


0.36 


< 9.6 x 10~ 5 


1.1 x 10~ 3 


TWA 6 


< 110 


> 2.6 


0.02 


0.14 


< 1.4 x 10~ 4 


1.4 x 10" 3 


TWA 7 


80 


6.8 


2. 1 


0.31 


2.0 x 10 -3 


7.8 x 10 -4 


TWA 8A 


< 150 


> 1.5 


0.04 


0.19 


< 6.6 x 10~ 4 


1.3 x 10 -3 


TWA 8B 


< 70 


> 3.2 


0.19 


0.04 


< 7.2 x 10~ 4 




TWA 9A 


< 60 


> 8.6 


1.3 


0.15 


< 6.8 x 10~ 4 


1.3 x 10~ 3 


TWA 9B 


< 60 


> 4.2 


1.5 


0.03 


< 3.2 x 10~ 3 




TWA 10 


< 80 


> 3.7 


0.09 


0.09 


< 2.6 x 10~ 4 


1.1 x 10~ 3 


HR 4796A 


108 


30 


110 


19.7 


4.8 x 10~ 3 




TWA 12 


< 65 


> 6.0 


0.94 


0.09 


< 1.0 x 10~ 3 


1.0 x 10 -3 


TWA 13A 


65 


7.8 


1.1 


0.18 


8.6 x 10" 4 


1.8 x 10" 3 


1 WA lO-D 


DO 


I . f 


"1 A 
1.1 


U.l I 


o.y X 1U 




TWA 14 


< 65 


> 5.1 


0.59 


0.07 


< 8.8 x 10~ 4 


1.4 x 10~ 3 


TWA 15A 


< 60 


> 3.5 


0.27 


0.03 


< 8.2 x 10~ 4 


4.7 x 10" 3 


TWA 15B 


< 65 


> 3.1 


0.20 


0.03 


< 7.9 x 10" 4 




TWA 16 


< 80 


> 4.5 


0.15 


0.12 


< 2.8 x 10~ 4 


7.0 x 10" 4 


TWA 17 


< 60 


> 5.1 


0.54 


0.06 


< 7.8 x 10~ 4 


9.8 x 10~ 4 


TWA 18 


< 80 


> 3.1 


0.06 


0.06 


< 2.6 x 10~ 4 


9.8 x 10" 4 


TWA 19A 


< 95 


> 13 


1.0 


2.10 


< 2.2 x 10~ 4 


7.1 x 10~ 4 


TWA 19B 


< 70 


> 11 


2.7 


0.38 


< 8.2 x 10~ 4 





a The minimum distance from a star for dust in thermal equilibrium at temperature, Tpds, is given by Rpds = 
0.5(T* /Tpos) 2 R*, where R* is radius of the star. 

b The estimated minimum mass in dust, Md, consistent with the observations for systems where Tpds can be determined or 
constrained, and assuming that the circumstellar debris is optically thin. A n average dust grain size of 2.8 /im in radius and a 
density for the material of p = 2.5 g cm" 3 are assumed iChen fc Jural 1200 If) . For systems with no detected IR excesses, Md is 
an estimate using the limits determined for Rpds and Lip/L,. 

C A distance to the TWA of 55 pc is assumed except for the five systems with measured parallaxes. For these systems, the 
distances listed in Table 1 are used to determine L* . 

d TWA 2, Hen 3-600, and HD 98800B are binaries and both stellar components are included in the luminosity values. 
G The minimum dust mass for the HD 98800B debris system is likely to be underestimated bec ause the assumption that the 
circumstellar material is optically thin in the IR is not valid fsee. lLow. Hines. fc Schneiderlll99fl(l . 
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Fig. 1. — MIPS 24 (im photometry compared to the 2MASS _ft" s -band measurements of the TWA sample. The dashed line represents the 
^24nm-to-i ? x s ratio under the assumption that both photometric bands lie on the Raylcigh-Jeans tail of the stellar spectrum. 
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Fig. 2. — The ratio of the observed 24 fim flux density to that determined for the _K" s -band from the 2MASS Point Source Catalog plotted 
against the best-fitting model temperature of the stars. TW Hya, Hen 3-600, HD 98800, HR 4796A, and TWA 19A are not shown so that 
TWA members with T < 6000-ft" or without a large 24 (im excess can be displayed in more detail. The solid curve represents the flux ratio 
between the two filter bandpasses for a blackbody. Likewise, the dashed curve is the ratio determined from the NextGen models used to 
estimate the 24 (im and 70 (im photosphcric flux densities. 
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Fig. 3. — The 70 /im excesses measured by Spitzer relative to the predicted 70 fim flux densities of the stars in the TWA. As in Figure 2, 
the data are plotted against the model temperature of the stellar photospheres. Each star detected by MIPS at 70 fim is labeled. One-sigma 
upper limits are denoted by open triangles. 
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Fig. 4. — The SEDs of the six TWA members with confirmed IR excesses. Model stellar photospheres are shown for the stars, along 
with available photometric measurements from 0.36-160 fim. Measurements made by MIPS are denoted by filled circles. Data from the 
literature are shown by open symbols. References to the optic al (circles), near-I R (squares), ground-based mid-IR (triangles) , and IRAS 
(diamonds and upper limits) photometry are given in Table 1. lUehrz et aT] (1999) showed that HD 98800B produces the IR excess in this 
binary system, though the sum of the models and observed flux densities of A and B are plotted. Blackbody fits to the MIPS and IRAS 
data for HD 98800B and HR 4796A are shown as bold dashed curves. The bold solid curves shown for TWA 7 and 13AB are the sum of 
the photospheres and IR excesses assuming that the excesses are described by blackbodies having Tpris = 60—80 K. Color temperatures 
(T c ) for TWA 7 and 13AB are also displayed. The SEDs have been spaced apart for clarity and a 1000 mjy fiducial is labeled for each 
star. To recover the measured flux densities in mjy, multiply the SEDs by the following factors: 0.33 (TW Hya), 0.29 (Hen 3-600), 0.67 
(HD 98800AB), and 0.2 (HR 4796A, TWA 7, and 13AB). 
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Fig. 5. — MIPS SED-modc spectra of TW Hya and HD 98800B. For comparison, the spectrum of a Boo is displayed and was used 
to calibrate the fluxes and spectral shapes of the two TWA members. The broad-band flux densities measured from the MIPS 70 fim 
photometry are shown for TW Hya (square) and HD 98800B (triangle). The spectra are calibrated assuming that F v = 14.7 Jy at the 
effective wavelength of the 70 fim bandpass for a Boo (circle), and that its spectrum is adequately described at these wavelengths by 
The best-fit blackbody to the IR excess of HD 98800B using the broad-band infrared measurements from IRAS, Spitzer-MIPS, 



and the ground is shown by the dashed curve and has a temperature, Tbb = 160 K (see Figure 4). 
SED-mode spectra is caused by pixel-to-pixel sensitivity variations of the 70 /im array. 



Most of the apparent "noise" in the 



